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X(5) is a paradigm for the structure at the critical point of a particular first-order 
phase transition for which the intrinsic energy surface has two degenerate minima 
separated by a low barrier. For a finite system, we show that the dynamics at such a 
critical point can be described by an effective deformation determined by minimizing 
the energy surface after projection onto angular momentum zero, and combined with 
two-level mixing. Wave functions of a particular analytic form are used to derive 
estimates for energies and quadrupole rates at the critical point. 
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Quantum phase transitions, occurring at zero temperature as a function of a coupling 
constant, have become a topic of great interest in different branches of physics [1, 2, 3]. 
Often this type of phase transition involves a structural change between different shapes or 
geometric configurations. Advanced experiments have identified such quantum shape-phase 
transitions in a variety of mesoscopic systems, with a finite number of constituents, e.g., 
nuclei, molecules and atomic clusters. A key issue is to understand the modifications brought 
in by the finiteness of these systems near criticality. This question can be conveniently 
addressed in a class of quantum models in which the Hamiltonian is expanded in elements 
of a Lie algebra. Such models are widely used in the description of nuclei and molecules [4, 5]. 
In the present work we study this question in connection with nuclei, exemplifying a finite 
system undergoing a first-order shape-phase transition. The tools developed are applicable 
to other mesoscopic systems described by similar models albeit with different spectrum 
generating algebras. 

Recently, it has been recognized that quantum shape-phase transitions are amenable to 
analytic descriptions at the critical points [1, 2]. The importance of these analytic bench- 
marks of criticality lies in the fact that they provide a classification of states and analytic 
expressions for observables in regions where the structure changes most rapidly. For nu- 
clei these benchmarks were obtained in the geometric framework of a Bohr Hamiltonian for 
macroscopic quadrupole shapes. In particular, the E(5) [1] (X(5) [2]) benchmark is applica- 
ble to a second- (first-) order phase transition between spherical and deformed 7-unstable 
(axially-symmetric) nuclei. In the present work we focus on the X(5) benchmark for which 
an empirical example has been found in 152 Sm, 150 Nd, 156 Dy and 154 Gd [6, 7, 8, 9]. The 
role of a finite number of nucleons can be addressed in the algebraic framework of the in- 
teracting boson model (IBM) [4] which describes low-lying quadrupole collective states in 
nuclei in terms of a system of N monopole (s) and quadrupole (d) bosons representing va- 
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lence nucleoli pairs. The model has U(6) as a spectrum generating algebra and its three 
dynamical symmetry limits: U(5), SU(3), and 0(6), describe the dynamics of stable nuclear 
shapes: spherical, axially-deformed, and 7-unstable deformed. Phase transitions for finite 
N are studied by an IBM Hamiltonian involving terms from different dynamical symmetry 
chains [10]. An important conclusion from many such studies [11, 12, 13, 14, 15, 16, 17] is 
that although in finite systems the discontinuities at the critical point are smoothed out, 
features of the phase transition persist even at moderate values of iV (N = 10 for 152 Sm). 

The original formulation of the E(5) and X(5) models [1, 2] employed an infinite square- 
well potential in the (5 variable of the Bohr Hamiltonian. This is an adequate approximation 
for an E(5)-like second-order phase transition, where the relevant potential is 7-indpendent 
and flat-bottomed. This behaviour persists in the finite-N energy surface, obtained by 
the method of coherent states [10, 18]. In this case, the large fluctuations in (3 can be 
taken into account by means of an effective /3-deformation determined by minimizing the 
energy surface after projection onto the appropriate [0(5)] symmetry [14]. The structure 
of X(5)-like first-order phase transition is more complex. The geometric X(5) solution [2] 
assumes the decoupling of the (3 and 7 degrees of freedom. The corresponding finite-N energy 
surface displays two degenerate minima separated by a low barrier. Analysis of the IBM 
wave functions near the critical-point show clear evidence for phase coexistence and level 
crossing [11, 13, 15]. This indicates that although the barrier between the two minima is 
small, its effect cannot be completely ignored. In the present work we examine the conditions 
for and properties of A(5) critical-point structure in a finite system. 

Focusing on the dynamics of the (3 degree of freedom, in the geometric approach, the X(5) 
eigenfunctions [2] are proportional to Bessel functions, and the spectrum consists of families 
of states, Lig, labeled by £ = 1, 2, . . ., with angular momentum L — 0, 2, 4, . . . and projection 
K = along the symmetry-axis. The X(5) benchmark leads to analytic parameter-free 
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predictions for energy ratios and B(E2) ratios which, as seen in Table I, are in-between the 
values expected of a spherical vibrator [Z/ (5)] and an axially-deformed rotor [577(3)]. 
In the algebraic approach, the £7 (5)-S77(3) transition is modeled by the Hamiltonian 

H = eh d -nQ-Q . (1) 

Here h d is the <i-boson number operator, Q is the quadrupole generator of S77(3) and the 
dot implies a scalar product. In the U(5) limit (k = 0), the spectrum of H is harmonic, 
and the eigenstates are classified according to the chain £7(6) D £7(5) D 0(5) D 0(3) 
with quantum numbers | N, rid, t, L) (for r > 6 an additional multiplicity index may be 
required for complete classification). In the S77(3) limit (e = 0), the Hamiltonian is related 
to the Casimir operator of SU (3) and the spectrum is indicated in the caption of Table I. 
The eigenstates are classified according to the chain £7(6) D S77(3) D 0(3) with quantum 
numbers | N, (A, /i), K, L). A geometric visualization is obtained by an intrinsic energy 
surface defined by the expectation value of the Hamiltonian in the coherent state [10, 18] 

\(3, r ,N) = (An)-V2 (6 tyv| ), (2) 

where b\ = (1 + (3 2 )-^ 2 [(3 cos-f d\ + (3 sin 7 (4 + dL 2 )/y/2 + st] with (3 > and < 7 < tt/3. 

The IBM Hamiltonian at the critical point of the U (5)-577 (3) phase transition corresponds 
to the following choice of parameters [10] in the Hamiltonian of Eq. (1) 

c = £«(2JV-3). (3) 

Its intrinsic energy surface given by 

E(P, 7 ) = -5nN + KN 2 ^~^f (l " 4^/3 cos 3 7 + 8/3 2 ) , (4) 

has the typical form of a Landau potential for a first-order phase transition, with two 
degenerate minima, at f3 — and at ((3 = 7775,7 = 0). As shown in Fig. (1), the barrier 



separating the spherical and prolate-deformed minima is extremely small and the resulting 
surface, E{(3) = E(f3, 7 = 0), is rather flat. This behaviour motivated the use of a square-well 
potential in the X(5) model [2]. Guided by the experience gained with finite-N flat-bottomed 
potentials in second-order phase transitions [14], we are led to consider states, \(3;N,L), of 
good 0(3) symmetry L projected from the intrinsic state 7 = 0; N) of Eq. (2), with an 
effective /5-deformation yet to be determined. In the U (5) basis these L-projected states are 

I P\ N,L) = £ \ [l + (-1)—] e d f\ N, n d , r, L> , (5) 

where L is even, r takes the values compatible with the 0(5) D 0(3) reduction and the n d 
summation covers the range r < n d < N. The coefficients are of the form 



^ nd ' T L N Jt If AT \\f \\\f , , Q\||ll/2 ^ > 

[(N — n d )\[n d — T)\\(n d + r + 3)!!J 
where ^((3) is a normalization factor. In some cases analytic expressions of these coeffi- 
cients can be derived. Specifically, for L = 0: /<°) = (-l) T v / 2r + 3 with r = 0, 3, 6, . . . 
and for L = 2: f& = (= (-1) T+1 V7+T) with r = 1,4,7,... (with 

r = 2, 5, 8, . . .). In general, can be obtained by numerical diagonalization of Hamil- 

tonians which have the condensate (2) as an exact eigenstate [19]. The states \(3;N,L) 
interpolate between the U(5) spherical ground state, \s N ), with n d = r = L = 0, at (3 = 0, 
and the 577(3) deformed ground band with (A, //) = (27V, 0), at f3 — y/2. For arbitrary (3 the 
matrix elements of the Hamiltonian in these states define an L-projected energy surface 
E^\f3) = ((3; N, L\H\(3; N, L), which can be evaluated in closed form 
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The quantities S[ N ^ , S^l and E^l are, respectively, the expectation values of n s = s^s, 
(s t ) 2 s 2 and n s h d in the states \(3;N,L). They are determined by the normalization factor 



r%\p) ofEq. (6) 
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with Sg» = SS'SS"" and E<5 = (JV - l)Sg> - S™. 

As noted in [21] and shown in Fig. (lb), the L = projected energy surface, E^} ((3), no- 
longer exhibits the double minima structure observed in the (unprojected) intrinsic energy 
surface. Instead, there is a minimum at (3 > 0, a maximum at (3 — 0, and a saddle point in 
the 7 direction at (3 < 0. E^ 2 (f3) resembles the potential used in the geometric collective 
model calculation of 152 Sm [22] with a minimum at a larger value of (3 > 0, and a flat shoulder 
near (3 = 0. The different behaviour of E^} 2 ((3) and -E1=o(/3) can be attributed to the fact 
that, as shown in Fig. (la), the L = 2 state is well above the barrier and hence experiences 
essentially a flat-bottomed potential. In contrast, the two minima in the intrinsic energy 
surface support two coexisting L = states which, in view of the low barrier, are subject to 
considerable mixing. The mixing between the spherical and deformed L = states 

I0i) = \s N ) , 

\<h) = \(3;N,L = 0) , (9) 
can be studied by means of a 2 x 2 potential energy matrix, = {(f>i\H\(pj) , given by 
m n = -5k N , mi2 = -kN \(3 2 (N - 1) + 5 



r\2 , 



m 22 = e { l n Jm , r 12 = [mr^'mi " , (io) 

where r 12 = (0i|02) is the overlap. The eigenpotentials are obtained by diagonalization and 
the equilibrium /3-deformation is determined from the minimum of the lowest eigenpotential. 
This procedure is in the spirit of the matrix coherent-state approach proposed in [23] to 
describe the geometry of configuration mixing in nuclei. An important difference being that 
instead of using coherent states with different values of N, in the present study we consider 
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states with good angular momentum, L = 0, and fixed N. A slight complication arises from 
the fact the two states |0i) and \(f) 2 ) are not orthogonal. This can be taken into account by 
transforming into an orthonormal basis of L = states 

l*i> = l<M , 

l* 2 ) = (l-rl^d^-nM) . (li) 

The matrix elements = (^i\H\^j) are then given by 

Ku = mn , K12 = (1 - rf 2 y 1/2 ( "I12 - r 12 m u ) 

K 22 = (1 - r\ 2 y l (m 22 - 2r 12 m 12 + r\ 2 m X \ ) . (12) 
The eigenvalues define the eigenpotentials 

E^lM = (K 11 + K 22 ±A)/2, 



A = ^{K 22 -K u y + AKt 2 , (13) 

and the corresponding eigenvectors are given by 

|$£2o) = sin 9 l^i) + cos 9 |* 2 ) , 
|$« ) = cog^i^) - S in^|^ 2 ) , 

The lowest eigenpotential, E^} Q (P), is shown in Fig. (lb) for N = 10. It has a minimum at a 
value of (5 — 0.591, larger than the minimum (f3 = 0.528) of the (unmixed) L = projected 
energy surface and in close proximity to the minimum (f3 = 0.619) of the L = 2 projected 



energy surface. We now identify the members of the ground band (£ = 1) as (0^) = |$ 
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L=0/> 



Eq. (14), with energy E ( L J (p), Eq. (13), and for L > even, |L+) = N, L), Eq. (5), with 
energy E^\(3), Eq. (7). The bandhead of the first excited band (£ = 2) corresponds to 



|0+) = l^o), Eq. (14), with energy E ( L + J {l3), Eq. (13). The value of f3 used in the indicated 
wave functions and energies is chosen at the global minimum of E^J (P). 

Having at hand explicit expressions for the wave functions, we can evaluate E2 matrix 
elements. To conform with the geometric X(5) model [2], we employ an IBM quadrupole op- 
erator, T(E2) = Ss+s^d, linear in the deformation. For transitions involving the L = Of, 2f 
and Of states, analytic expressions can be derived for B(E2) values by means of the matrix 
elements T x = (f3;N,L' = 2\\T(E2) \ \(3; N,L — 0) and T 2 = ((3; N, L' = 2\\T(E2) \\s N ), 

j-i — - 



T2 = 7 i%72 • ( 15 ) 

To test the suggested procedure we compare in Table II the U(5) decomposition of exact 
eigenstates obtained from numerical diagonalization of the critical Hamiltonian, Eq. (3), 
for N = 10, with that calculated using the L-projected states with (5 = 0.591 [the global 
minimum of E^ L J (/3)]. As can be seen, the latter provide an accurate approximation to 
the exact eigenstates for yrast states (the overlaps between calculated and exact states 
are 99.6, 99.1, 98.0, 97.0, 96.6, 96.65% for L = Of, 2f, 4f, 8f, lOf, respectively). This 
agreement in the structure of wave functions is translated also into an agreement in energies 
and B(E2) values as shown in Table I. For the L = Of state, the agreement is fair but 
less precise, in view of the smaller overlap (81.9%). It appears that this non- yrast state is 
affected by additional states beyond the two-state mixing considered here. The results of 
Table I and II clearly demonstrate the ability of the suggested procedure to provide faithful 
estimates to the exact finite-N calculations of the critical IBM Hamiltonian, which in-turn 
capture the essential features of the X(5) critical-point structure relevant to 152 Sm (which 
is slightly past the phase transition towards SU(3) [6]). 

To summarize, we have examined the structure at the X(5) critical-point of a first- 
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order shape-phase transition in a finite system, by means of an effective /3-deformation, 
determined by variation after L-projection and two-level mixing. The same procedure can 
be used throughout the coexistence region, where the two minima in the potential coexist 
but are not necessarily degenerate. In this case, the control parameter n/e in Eq. (1) 
spans the range between the spinodal point (where the second minimum shows up) and the 
antispinodal point (where the first minimum disappears), which embrace the critical point. 
Although we have treated explicitly the case of atomic nuclei described by the IBM, the 
same approach is applicable to first-order phase transitions in other mesoscopic systems, 
such as polyatomic molecules described by the algebraic vibron model [5]. This work was 
supported by the Israel Science Foundation. 
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TABLE I: Excitation energies (in units of E(2f) = 1) and B(E2) values (in units of B(E2; 2+ -> 
0+) = 1) for the X(5) critical-point benchmark [2], for several N=10 calculations, and for the 
experimental data of 152 Sm [6]. The finite-N calculations involve the exact diagonalization of the 
critical IBM Hamiltonian [Eq. (3)], 7-projected states [Eqs. (7), (13) and (15) with (3 = 0.591], the 
17(5) limit [en d ] and the S77(3) limit (k/2)[-A 2 - fi 2 - Xji - 3A - 3/i + 37(L + l)/4]. 
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18.33 


13.21 




5.67 


2.64 
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25.33 


5.62 
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1.58 


1.61 
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1.8 
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1.98 


1.85 
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1.48 
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2.27 


1.92 
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2.61 
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1.37 


2.22 
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0.78 


0.61 
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TABLE II: U (5) decomposition (in %) of yrast Lf states for iV = 10. The exact values are obtained 
from numerical diagonalization of the critical IBM Hamiltonian, Eq. (3). The calculated values are 
obtained from the L-projected states, Eqs. (6) and (14) with (5 = 0.591. 
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FIG. 1: Energy surfaces of the critical IBM Hamiltonian, Eq. (3), with k = 0.2 and ./V = 10. 
(a) Intrinsic energy surface E(j3) = E((3,j = 0), Eq. (4), [solid line]. The unmixed L = and 
L = 2 levels are shown for illustration, (b) E{(5) [solid line] as in (a), unmixed L = [dashed line] 
and L = 2 [long dashed line] projected energy surfaces, £x(/3) = E^\f3), Eq. (7), and the lowest 
L = eigenpotential [dot-dashed line], E^2 Q (f3), Eq. (13), whose global minimum is at f3 = 0.591. 
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